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Abstract Recently, S.W. Kahler studied the solar energetic particle (SEP) event timescales 
associated with coronal mass ejections (CMEs) from spacecraft data analysis. They obtained 
different timescales of SEP events, such as TO, the onset time from CME launch to SEP 
onset, TR, the rise time from onset to half the peak intensity (0.5/p), and TD, the duration 
of the SEP intensity above 0.5/p. In this work, we solve SEPs transport equation considering 
ICME shocks as energetic particle sources. With our modeling assumptions, our simulations 
show similar results to Kahler’s spacecraft data analysis that the weighted average of TD 
increases with both CME speed and width. Besides, from our simulation results, we suggest 
TD is directly dependent on CME speed, but not dependent on CME width, which were not 
achieved from the observation data analysis. 
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1 INTRODUCTION 


Solar energetic particle (SEP) events could be mainly divided into two classes through duration and in¬ 
tensity. The short-duration and low-intensity events, which are called impulsive events, are considered to 
be produced by solar flares. On the other hand, the longer duration and higher intensity ones, which are 
called gradual events, are considered produced by coronal and interplanetary shocks driven by coronal 
mass ejections (CMEs). It is interesting to study the relationship between gradual SEP event pr operties and 


the characteristics of the associated CMEs. With the first-order Eermi acceleration mechanism (iZank et al 


2000h introduced an onion shell model using a one-dimensional hydrodynamic code for the evolution of 


the CME-driven shock in the Parker interplanetary magnetic field (IM F). The model is valid only in strong 


shocks due to Bohm diffusion coeffi cient used, so 
strengths. In addition. 


Li et al. 


Rice et al 


(l2003h modified it to be usable in arbitrary 


(120031) studied the transport of SEPs with their onion shell acceleration model 


considering particles pitch angle scattering without perpendicular diffusion. In their model, charged parti¬ 
cles’ pitch angle diffusion is not considered between two consecutive pitch angle scatterings. Furthermore, 
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Verkhoglvadova et alJ (120091120101) adopted this model to study individual SEP events caused by CME 
shocks, their simulation results can fit well with spacecraft observations for different elements. On the other 
hand, considering that the interpla netary coronal mass ejection (ICME) shocks can cont inuously accelerate 


SEPs when propagating outward. 


Kallenrode & Wibberena (119971) : 


shock as a moving particle source. And the model was adopted in a numerical codePIbv iWang et al 


Kallenrode ( 2001 ) treated the ICME 

T 


( 120121 ) 


to study ICME driven shock accelerated particles’ transport in three dimensional solar wind and IME includ¬ 
ing both parallel and perpendicular di ffusion co e fficien ts. Eurthermore, under varying perpendicular diffu¬ 
sion and shock acceleration strength. 


Pin et al 


(120 131) reproduced the reservoir phenome non with SPTC 


numerical simulations. In addition, with the same numerical modeling. 


Wang & OinI (12015) researched the 


gradual SEP events spectra forcusing on the spatial and temporal invariance. Einally, 


Pin & Wangl (12015h 


compared the simulation results from SPTC with the multi-spacecraft {Helios 1, Helios 2, and IMP 8) 
observations during a gradual SEP event, and they obtained the SPTC simulations which best fit the SEP 
event observed by spacecraft located in different space. 

To investigate the relationship between SEP event properties with the associated CMEs, 


Ding et al. 


(120141) studied the interaction of two CMEs erupted nearby during a large SEP event by multiple spacecraft 
observations with the graduated cylindrical shell model. And they obtained the solar particle release time 


and p ath length which indicated the necessary influence of the ”twin-CME” (ILi et al 


2012 ; 


Temmer et al 


2012) on the SEP event. 


Bec ause of the huge d amage caused by SEPs, the study of peak intensities of SEPs becomes very im¬ 


portant. 


Ding et al 


(12015h presented the new observation r esults of peak intensi t y with Ee/P ratio, which 


indicate the role of seed population in extremely large SEPs. 


Reinard & AndrewsI (120061) studied the depen 


dence of the occurrence and peak intensities of SEP events with CME properties thoroughly using databases 
of the LASCC/5 OHO CMEs and the GOES E > 10 MeV protons. Besides peak intensities, timescales are 
another very important property of SEPs which could make contribution to both space weather forecasting 
and understanding of the SEP injection profiles and propagation characteristics. 

In order to study the properties and associations of SEP events, [j 


Cane et al 


(1201 Oh compared SEPs with 


flares and CMEs of 280 solar proton events which extended above 25 MeV occurred from 1997 to 2006 by 
near-Earth spacecraft. They divided the events into 5 groups according to the ratios e/p and Ee/O at event 
onset. Their results suggested that SEP event occurrence and peak intensities are more likel y to be associated 


with faster and wider CMEs, especially with western CME source regions. Eurthermore, 


Pan et al 


(1201 ih 


investigated SEP times cales, such as the SEP onset time, the SEP rise time, and the SEP duration. With an 
ice-cream cone model. 


Pan et al. 


(1201 Ih studied LASCO/5 OHO observation data of 95 CMEs associated 


with SEP events during 1998 - 2002, and came to conclusions that the SEP onset time has no significant 
correlation with the CME speed, nor with the CME width. They also suggested that the SEP rise time and 
the SEP duration have significantly positive correlations with the radial speed and angular width of the 
associated CMEs unless the events are not magnetically well connected to the Earth. 


Kahlen (120 13h did a research on the relationship between the EPACT/Wind 20 M eV SEP event s 


timescales and their associated CME speed and widths observed by LASCO/5 OHO. In 


Kahled (l2013h . 


* Hereafter, we denote tFie code as Slioclc Particle Transport Code, SPTC. 
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217 SEP events observed in a solar cycle during the period 1996-2008 were used. They defined the three 
characteristic times of the SEP events. The time from inferred CME launch at 1 Rq to the time of the 20 
MeV SEP onset at Wind was denoted as TO. The time from SEP onset to the time the intensity reached half 
of the peak value (0.5/p)was denoted as TR. And the time during which the intensity was above Q.SIp was 
denoted as TD. Prom their results, they found that CME speed and width were of signihcant correlation and 
it is not easy to interpret the contribution of CME speed and width to timescales separately. Therefore, they 
suggested that faster and wider CMEs which drive shocks and accelerate SEPs over longer times would 

thus produce the longer SEP timescales TR and TD. _ 

In this paper, with the data used in the analysis of iKahlen (120 13h . we stud y the CME tim escales by 


numerical simulations with the SPTC, and we compare our results with that of 


Kahlerl (120 13h . In section 


2, we present the model. In section 3, we present the data analysis. In section 4, we show our results. In 
section 5, we present the conclusions and discussion. 


2 MODEL 


We m odel the transport of SEPs by following previous research ( e.g.. 


2009 ). The three-dimensional focu sed transport equation is written as (ISkillin 


Pin et al 


2006; 


Zhang et al 


20091) 


Qinet^ 


3 


1971 


2006; 

Zhan 

et al. 

Schlickeisei 

2002 
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where f{x,p, p, t) is the gyrophase-averaged distribution function, x is the position in a non-rotating helio¬ 
graphic coordinate system, p is the particle pitch-angle cosine, p is the particle momentum, v is the particle 
speed, t is the time, and are the particle perpendicular and pitch-angle diffusion coefficients, re- 

A /A 

spectively, y™ = y™r is the solar wind velocity which is in the radial direction, and L - (ft VlnBo) 
is the magnetic focusing length determined by the magnitude of the background magnetic field Bq and the 

A 

unit vector along the local magnetic held ft. In the equation ([1]), almost all important transport effects are 
included, i.e., perpendicular diffusion (1st term in RHS), pitch angle diffusion (2nd term in RHS), particle 
streaming along held line and solar wind bowing in the IMP (third term in RHS), adiabatic cooling in the 
expanding solar wind (4th term in RHS), and magnetic focusing in the diverging IMP (5th term in RHS). 
Here, the drift effects are neglected for lower-energy SEP transport in the inner heliosphere. Also the IMP 
is modeled with the Parker held. 

By following [ 


Burger et al. 


(120081), diffusion coefficients are dete rmined. We set the per pendicular dif¬ 


fusion coefficient from the nonlinear guiding cente r (NLGC) the ory dMatthaeus et al 


with the analytical form according to 


Shalchi et al 


(12004 


2010I) . 


20031) approximated 


„ _„/2/3A/3 

A-j. - Vl. A.. 


(i-bb). 


( 2 ) 


where is a parameter to control the value of the perpendicular diffusion coefficient. Por simplicity, 
Aj^ is set to be independent of p with the assumption that particle pitch-angle diffusion is much faster than 
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perpendicular diffusion , but generally dependent perpendicular diffusion coefficient should be used (e.g.. 


Pin & Shalchi, 


20141 


3 . 


The parallel 

1968 

Earl 

1974) 


particle mean free path (mfp) /l|| is written as (iJokiDii . 


1966, 


Hasselmann & Wibberenz 


^11 


T J-i 


2\2 


D„ 




(3) 


-1 infill 

and parallel di ffusion coefficient /tn can be writ ten as /ch - V/lii/ 3 . 


We follow 


Beeck & Wibberena(ll986ll and 


Teufel & Schlickeiserl (120031) to model the pitch angle diffu¬ 


sion coefficient 

{|;U I '-' + h} (1 - , (4) 

where G is a parameter to control the value of v is the particle speed, Ri = pc!(\q \ Bq) is the particle 

Larmor radius. Here, a larger value of h - 0 . 01 is c hosen for non-linear effect of pitch angle diffusion at 


ju = 0 in the solar wind (lOin & Shalchil 


2009 


201 4ah. 


To model the particle injectio n, the shock is treated as a moving SEP source with the boundary condition 


dKallenrode & Wibberenz , 


1997): 


/ r 

10(6', <f>)\ 

1 -) 

4>cip) 


p - \(l)(,6,ip)\) 


(5) 


a - < 


where a and (fee are the shock acceleration strength parameters. We assume (fee as a constant, but a as a 
function of shock speed, e.g., we set 

-3.5 if Vs < vi 

Vj/vo - 5 if vi < Vs < V2 (6) 

^ -2 if Vs > V2 

where Vo = 500kms“',vi = 750 km s“', and V 2 = 1500 km s“'. 0(0, is the angle between source center 
and any point of particle injection (0, ip). y is the spectral index of source particles. In the simulations, we 
inject energetic particle shells with small space intervals Ar. H{x) is the Heaviside step function, with 0j 
being the half ang ular width o f the sh ock. A more detailed description of the shock model of our simulations 


can be referred to 


Wang et af 


mm . 


The transport equation ([TJ is solved by a time-backward Markov stochastic pro cess method (IZhang , 


1999h in the simulations. And the detailed description of the method can be referred to 


Pin et af 


(l2006h . As 


mentioned in section [T] our numerical code of transport of energetic particles with the CME driven shock 
as a moving particle source is denoted as Shock Particle Transport Code, i.e., SPTC. 

3 DATA ANALYSIS 

We investigate 20 MeV proton intensity-time profiles of SEP events during 1996 to 2008 with their associ¬ 
ated CMEs. In particular, the SEP data is from EPACT (the Energetic Particles: Acceleration, Composition, 


and Transport) (Ivon Rosenvinge et af 


1995h experiment on the Wind spacecraft, and the information of 


their related CMEs is observed by S OHO (th e S olar and Hel i osphe ric Observatory mission) LASCP 


(Large Angle and Sp ectrometric Coronagraph) dBrueckner et al 


period (iKahlei 


19951) . Pf the total 217 SEPs during this 


20131) . we study 204 SEPs whose CME parameters are available. In addition, for each event. 


the CME sola r source is dete rmined by flare location, and the speed {vcme) and width {Wcme) of CME are 
obtained from Kahler ( 2013 ). 
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3.1 Parameter Selection 


For the grouping and selecting data, we follow the method suggested by 


Kahleri (12013h dividing the 204 


events into five longitude ranges with about 41 events each, and subdividing each longitude range into 
several groups sorted on vcme and Wcme, respectively. The median values of longitude, vcme, and Wshock 
in each group are used as the characteristic values. 

From data analysis of spacecraft observations, it is not easy to identify SEP onset time accurately which 
is usually covered by the background of intensity. Therefore, in this work, we only focus on the variation 
of TD with Vcme and Wcme- To compa re with the obs ervation, we obtain the data analysis results of vari- 

Kahlen (120 13h as shown in Table [T] From Table [T] we can see, we 


ation of TD with vcme and Wcme from 


study SEP events with source location longitude in three ranges, W33-W60, W62-W90, and WlOO-bWL, 
with median values W48, W77, and W112, respectively. Note that bWL indicates sources behind the west 
limbs. In each range of longitude, TD is shown as varying with the median values of vcme and Wcme by 
subdividing the range into several groups sorted on vcme and Wcme, respectively. 

Table 1 The data analysis results of variation of TD with vcme and Wcme from lKahlen (120 13h . 


Source Location 

W33-W60 

W62-W90 

WlOO-bWL 

Longitude 








Vcme (km/s) 

TD (h) 

Vcme (km/s) 

TD (h) 

Vcme (km/s) 

TD (h) 

TD varying 

450 

6.3 

650 

6.5 

620 

13.2 

800 

12.0 

1150 

9.8 

900 

14.0 

with Vcme 

1175 

8.8 

1450 

21.3 

1325 

12.5 


1600 

14.5 

2100 

18.1 

1750 

17.0 


Wcme (°) 

TD (h) 

Wcme (°) 

TD (h) 

Wcme (°) 

TD (h) 

TD Varying 

77 

8.3 

133 

7.8 

100 

7.5 

with Wcme 

208 

11.3 

171 

15.0 

178 

13.3 


360 

15.8 

360 

16.4 

360 

17.2 


In order to study SEP timescales associated with CMEs, we use the SPTC described in Section |2] to 
simulate the transport of SEPs assuming the CME shock as a moving particle source and that the shock 
nose is in the flare direction relative to the solar center. In SP TC, the speed o f shock, v^, and the width of 
shock, Ws, are needed. While, in the spacecraft data analysis of lKahlen (120 13h the speed and width of CME 


are used instead. To compare the simulation results with the spacecraft data analysis, we need a model for 
relationship between and vcme, and that for relationship between Ws and Wcme- Eirstly, we assume the 
speed of CME is the same as that of shock, v., = vcme- Secondly, since the width of shock (W.,) is larger 
than that of CME (VTcM£),we set. 


Wcme + AIT if Wcme < 360° — AW 

W, = { (7) 

360° otherwise. 

By testing several value of AW, we finally set AW = 90°. It is noted that such kind of model for Ws is only 
an approximation, and it could lead to the discrepancy between the observation and simulation results. So 
we need to use a better W^ model in the future. Generally, the event source is near the solar equator, so 
the characteristic latitude of source location is set as 10° north. Other important simulation parameters not 
varying are shown in Table |2] 
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Table 2 Model Parameters Used in the Calculations. 


Parameter 

Physical meaning 

Value 

E 

Particles energy 

20 MeV 

ro 

Observer solar distance 

1 AU 

Ar 

Shock space interval between two fresh injections 

0.001 AU 

fc 

Radial normalization parameter 

0.05 AU 

7 

Spectral index of source particles 

-3.5 

</>c 

Shock strength parameter 

15° 


Particle mean free path 

0.16 AU 

Ke/kw"" 

Ratio between perpendicular and parallel diffusion coefficient 

6.1% 

fbO 

Inner boundary 

0.05 AU 

rti 

Outer boundary 

50 AU 


“ For 20 MeV protons in the ecliptic at 1 AU. 

In order to investigate the relationship between solar wind speed vsty and CME speed vcme, we obtain 
vsw observation data from Wind spacecraft for the 204 CME events to fit the relationship between vqme 
and vsvr- It is shown that vqme and vsw are positively correlated. As we assumed above that Vj = vcme, the 
relationship between vcme and vsw would turn to that between and vsw- Thus vsiv can be represented by 
Vj as 

vsvr = 1-77X lO^V+425, (8) 


here, vjvp and Vj are in the unit of km s'. We also divide the events into several groups sorted on vcme, and 
obtain the median values of vcme as the characteristic ones for each group. So we obtain the counterpart 
values Vj and vsw through the assumption above. And we use the characteristic ones in the simulations 
shown in Table [3] The Table |3] also shows the other input parameters in each sim ulation coming from the 
characteristic values of vcme, Wcme and source location longitude picked up fromy 
Tabled 


Kahlen (120131) shown in 


3.2 Simulation Output 

Eor each data point, 3200000 virtual particles are calculated in our simulations. In our simulations, we 
obtain the time profiles of SEPs with characteristic speed and width of CME, with which we can get the 
SEP timescale, TD. Eor example, in Eigure[T] we show simulation results of 20 MeV proton flux during an 
SEP event. In the simulation of Pigure[T] we set solar wind speed as 502.2 km s ', longitude as 48 degrees 
west, CME speed as 1600 km s ', CME width as 180°, other parameters are shown in Table|3 In EigurefT] 
the dotted line indicates the peak intensity {Ip) of the event, and the dash-dotted line indicates the half peak 
intensity. Ts and Tg indicate the earliest and latest time when the intensity is half peak, respectively. So we 
can obtain TD - Te - T^ from the time profile of intensity of simulation results. 

Erom the results of the simulations we can also get the weighted averages as following. Eor example, 
in each range of shock speed and longitude, we have three ranges of shock width, so we have three values 
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Table 3 Shock speed and width, and solar wind speed used in all simulations. 


Source location 

N10W48 

N10W77 

N10W115 

Vj (km/s) 

Vsw (km/s) 


Vj (km/s) 

Vsw (km/s) 

WAl 

Vj (km/s) 

Vsw (km/s) 

'T,C’) 

450 

433.5 

167 

650 

433.5 

223 

620 

433.5 

190 

450 

433.5 

298 

650 

433.5 

261 

620 

433.5 

268 

450 

433.5 

360 

650 

433.5 

360 

620 

433.5 

360 

800 

455.4 

167 

1150 

455.4 

223 

900 

455.4 

190 

800 

455.4 

298 

1150 

455.4 

261 

900 

455.4 

268 

800 

455.4 

360 

1150 

455.4 

360 

900 

455.4 

360 

1175 

444.1 

167 

1450 

444.1 

223 

1325 

444.1 

190 

1175 

444.1 

298 

1450 

444.1 

261 

1325 

444.1 

268 

1175 

444.1 

360 

1450 

444.1 

360 

1325 

444.1 

360 

1600 

502.2 

167 

2100 

502.2 

223 

1750 

502.2 

190 

1600 

502.2 

298 

2100 

502.2 

261 

1750 

502.2 

268 

1600 

502.2 

360 

2100 

502.2 

360 

1750 

502.2 

360 


20 MeV SEP events 



Fig.l Flux of 20 MeV proton during an SEP event with parameters shown in the text. The 
dotted line indicates the peak intensity of the event, and the dash-dotted line indicates the half 
peak intensity. and Tg indicate the earliest and latest time when the intensity is half peak, 
respectively. 


of TD from simulation results with same shock speed and longitude but different shock width. For the 
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three ranges of shock width we can get their percentage according to the number of events, with which the 
weighted value of YD is obtained from the individual values of TD. 

Further, w e study the relationship between CME speed and CME width using the observation data in 


Kahlen(l2013h . We subdivide each longitude range into several groups sorted on CME width. We get average 


CME speed for each group. The results are shown in Eigure|2]as the relationship between average of CME 
speed and the median value of CME width. The three data points in each longitude group of Eigurel^match 
those of the three CME width bins of Table [T] The line indicates fitting of the data. It is found that in 
statistics the average CME width increases with the increasing of CME speed. 


Correlations of speed with width 



CME WIDTH(DEG) 


Fig. 2 Variations of CME speed as a function of CME width with different source location. The 
crosses indecate CME speed averages of the CME width ranges in the source location range of 
W33-W60, and the squares are that of W62-W90, the triangles are that of WlO O-bWL. The li ne 


indicates fitting of the data. The symbols are from observation data analysed by 


Kahlej(l2013l) . 


4 RESULTS 


Figure [3] shows SEP timescale TD vs. CME speed for 20 MeV SEP events detected at 1 AU with different 
source locations in different pannels. The top, middle and bottom panels show different longitudes of source 
locations, 48° west, 77° west, and 115° west, respectively. The black squa res indicate spacecraft observation 
data in Table[T]which are obtained from the data analysis of lKahlen(l2013l) . The TD and CME speed for each 
data points correspond to those of Table [1] The red triangles indicate the weighted average of simulations 
according to the distribution of number of events with different CME widths for any given CME speed 
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interval obtained from the observation data in 


Kahletl (120 13b corresponding to the the abscissa of black 


squares. The red and black dashed lines indicate the linear fitting of the weighted average simulation results 
represented by the red triangles and that of the spacecraft observation data represented by the black squares, 
respectively. From Figure [3 we can see, the simulation results show the similar trend of observation data, 
that is, the SEP timescale TD increases with CME speed. 


Correlations of SEP Timescale TD with CME Speetd 
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Fig. 3 SEP event timescale TD vs. CME speed. Different pa nels indicate d ifferent source loca¬ 


tion. The black squares are from observation data analysed bv lKahleil (l2013h . The TD and CME 
speed for each data points correspond to those of Table [1] The red triangles indicate weighted 
average of simulation results. The black dashed lines indicate linear fitting of observation data. 
The red dashed lines indicate linear fitting of the weighted average of simulation results. 


Eigure|4]shows plot similar as Eigure[3]except that x-coordinate is CME width. The value for each black 
squares correspond to those of Table [T] The red triangles indicate weighted average of simulation results 
according to the distribution of number of events with dif ferent CME speeds for any given CME width 
interval obtained from the observation data in iKahleii (l2Q13h corresponding to the abscissa of black squares. 


Similarly as in Eigure[2 the red and black dashed lines indicate the linear fitting of the weighted average 
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of simulation results and that of the spacecraft observation data, respectively. From Figure |4] we can see, 
generally, the simulation results show the similar trend of observation data but with less slope, that is, the 
SEP timescale TD increases with CME width. However, from top panel of Eigure|4](N10W48) it is shown 
that, the observation results shows the SEP timescale TD increases with CME width, but the simulation 
results shows constant for different CME width. It is noted that our simulations could show deviation from 
observations due to modeling and statistical problems. 


Correlations of SEP Timescale TD with CME Width 
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Fig. 4 SEP event timescale TD vs. CME width. Different pa nels indicate d ifferent source loca¬ 
tion. The black squares are from observation data analysed bv iKahleii (120 13h . The TD and CME 


width for each data points correspond to those of Table [T] The red triangles indicate weighted 
average of simulation results. The black dashed lines indicate linear fitting of observation data. 
The red dashed lines indicate linear fitting of the weighted average of simulation results. 


Eor further study on the contribution of CME speed and width to timescales separately, we plot the 
individual simulations and weighted average of simulation results as follows. 
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Figure |5] shows simulations of SEP timescale TD vs. CME speed for 20 MeV SEP events detected at 
1 AU with different source locations in different pannels. Similar as Pigure[3j the top, middle and bottom 
panels show different longitudes of source locations, 48° west, 77° west, and 115° west, respectively. The 
yellow, green, and blue triangles indicate simulations with different CME widths corresponding to those of 
Table [1] The each data point of the yellow line shows a individual simulation with a distinct CME speed 
but a common CME width 133° and source location of N10W48, and so are the green and bule lines with 
other source location. Besides, the value of all data points are shown in Tabled The red triangles indicate 
the weighted average of simulations according to the distribution of number o f events with d ifferent CME 
widths for any given CME speed interval obtained from the observation data in iKahleii (120 13h . 



500 1000 1500 2000 

CME SPEED(KM/S) 


Fig. 5 Simulations of SEP event timescale TD vs. CME speed. Different panels indicate different 
source location. The yellow, green and blue triangles indicate simulations with different CME 
widths. The red triangles indicate weighted average of simulation results. 


Erom simulations of Eigure|5]we can see, every single colored line increases, that is to say, for the same 
CME width, TD generally increases with the increasing of CME speed, and the weighted average of TD 
from simulations also generally increases with the increasing of CME speed. Erom the other aspect, the 
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colored lines and symbols are almost overlap, from triangles with a common abscissa but different colors 
we can see, when CME speed is fixed, TD with different CME widths are almost same, meanwhile, when 
CME width is fixed, TD with different CME speeds are increased. So we su ggest fr o m our simulation that 


Kahlerl (120131) could not pick 


TD is dependent on CME speed but not on CME width, which analysis of 
out. 

Eigurel^shows plot similar as Eigure|5]except that x-coordinate is CME width. The yellow, green, light 
blue and purple triangles indicate simulations with different CME speeds corresponding to those of Table 
[T] The each data point of the yellow line shows a individual simulation with a distinct CME width but a 
common CME speed 450 km s ' and source location of N10W48 , and so are the green and bule lines with 
other source location. Besides, the value of all data points are shown in Table [3] The red triangles indicate 
the weighted average of simulations according to the distribution of number o f events with d ifferent CME 
speeds for any given CME width interval obtained from the observation data in lKahlen (120131) . 



100 150 200 250 300 350 

CME WIDTH(DEG) 


Fig. 6 Simulations of SEP event timescale TD vs. CME width. Different panels indicate different 
source location. The yellow, green, light blue, and purple triangles indicate simulations with 
different CME speed. The red triangles indicate weighted average of simulation results. 
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From simulations of Figure |6] we can see, the yellow, green, light blue and purple lines are almost 
aclinic, that is to say, for the same CME speed, TD generally keeps constant with different CME width. 
However, the red line which combines each individual line connecting data points of common CME speed 
simulations with weighted average increases, that is to say, for the same CME width, TD increases with the 
increasing of CME speed. In addition, the weighted average of TD increases with the increasing of CME 
width. The reason is that with larger CME width it is more likely that CME speed becomes larger, so the 
weighted average of TD becomes larger consequently. 


5 CONCLUSIONS AND DISCUSSION 


Generally, the accurate measurement of the first arriving particles in SEP events depends on the level of 
SEP flux background, so usually it is difficult to determine the timescales TR, TO, and TO + TR. However, 
TD, which indicates the duration of the SEP intensity above 0.5/p, has nothing to do with the first arriving 
particles, so the measurements of TD are relatively accurate. Therefore, we only study the timescale TD, 
but do not study TR, TO,orTO + TR. 

In this work, we use the SPTC to simulate the transport of SEPs assuming the ICME shock a s a moving 


particle source with parameters obtained from spacecraft observations analysed by 


Kahlerl (120 13h . and other 


parameters set as typical values of SEP events. Pr om simulation s we get SEP timescale TD and compare 
with TD valu es from spacecraft data analysis by 


Kahlerl (120131) . Prom spacecraft observations shown in 


Kahlen (120131) we obtain the contribution of CME speed with the same CME width, and we also obtain that 


of CME width with the same CME speed. Pinally, from simulation results of TZ) we can obtain the average 
of TD weighted with the observations contribution. 

Our simulations show that with the same CME speed, TD keeps constant with the increasing of CME 
width, but that the weighted a verage of TD increases with the increasing of CME width. Prom spacecraft 


data analysis in 


Kahlerl (l2013l) it is shown that TD, which is actually weighted average, increases with the 


increasing of CME width. In addition, our simulations show that with the same CME width, TD increases 
with the increas ing of CME sp eed, and that the average of TD increases with the increasing of CME width. It 


is also shown in 


Kahlen (1201 3h with spacecraft data analysis that the weighted average of TD increases with 


th e increasing of CME speed. Our simulations generally agree with spacecraft observations data analysis 


of 


Kahlerl ( 1201 3b that the weighted average of TD increase with both CME spee d and width. Pu rthermore, 


Kahlen (1201 3l) that TD is 


with our modeling assumptions, our simulations show some results not shown in 
dependent directly on CME speed, but independent on CME width. 

In order to study whether TD increases with CME width or speed by using observation data, one should 
choose SEP events with same CME speed but different CME width to show if TD increases with CME 
width, and also one should choos e SEP events with same CME width but different CME speed to show 
if TD increases with CME speed. 


Kahlen (120131) did not do it because of limitation of events number. But 


simulations do not have this limitation, and that offers us physical insigh ts behind the observations. We 


compare the weighted average of simulation to the resul t of 


our weighted average generally agrees with the result of 


Ka 


Kahlei 


ilei ( 20131) . and we can show the trend of 


l2013l) . so our work do not contradict the 
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observation result of 


Kahlerl (120 13h . Meanwhile, our individual results can be used to show if TD depends 


on CME width with same CME speed. 

The model we use to calculate flux includes many effects, such as the source, parallel and perpendicular 
diffusion, adiabatic cooling, etc., the overall effects could be very complicated, so we have to use numerical 
simulations to get the results. It is possible in some cases TD would decrease. But generally, TD has a trend 
to increase with the same CME width and increasing CME speed, and TD has a trend to be constant with 
the same CME speed and increasing CME width. Here, we compare the general trend between observations 
and simulations. 

We choose shock model conditions to favor larger particle injections with increasing speeds and widths 
in order to compare with observations. There are some parameters arbitrarily chosen and fixed in all sim¬ 
ulations, we tried different parameters, for example, we tested simulations with different value of shock 
strength parameter (pc, such as 10°, 15°, 18°, and 25°, and we found they would not change our general 
results. In the future, we would continue to study the parameter effects in our model. 

The observational evidences of the first detected SEP onse ts or relea ses associated with the good mag- 


netic connection to source were discussed in 


Ding et al 


(120161) . Besides, 


Rouillard et al 


(12011 


2012h sug¬ 


gested that SEP onsets co uld be considered as sociated with the modeled first connections of field lines to 
shocks. On the other hand. lOin & Wang|(l2015h showed the onsets from SPTC simulation results can fit well 


with that from observations of HELIOS 1, HELIOS 2, and IMP 8 at different longitudes simultaneously 
with perpendicular diffusion. It is interesting to compare the effects of these models carefully in the future. 
There are many aut hors working on numerical simula t ions to produce SEP profiles from the shock 


onion shell model (e.g.. 


Verkhoglvadova et al 


2009 


2010 ; 


Wang et al 


2012 


Pin et al 


20131) . and they 


usually study the individual SEPs in detail, in this work, however, we are trying to st udy many SEP s with 


simulations so we can compare with obseryations statistically. CME width data from 


Kahlen (120 13h were 


obseryed by only one satellite, S OHO, so they are lack of determinacy. In the future, we would study the 
CME data of multi-spacecraft obseryations. In addition, we would study peak intensity of gradual SEP 
eyents associated with CM Es by comparing the simulations of SPTC with the spacecraft data analysis (e.g.. 


Kahler & Vourlidas. 


20131). 
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